Band dispersion in the deep Is core level of graphene 
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Abstract 

Chemical bonding in molecules and solids arises from the overlap of valence electron wave func- 
tions, forming extended molecular orbitals and dispersing Bloch states, respectively. Core electrons 
with high binding energies, on the other hand, are localized to their respective atoms and their 
wave functions do not overlap significantly. Here we report the observation of band formation and 
considerable dispersion (up to 60 meV) in the Is core level of the carbon atoms forming graphene, 
despite the high C Is binding energy of ~ 284 eV. Due to a Young's double slit-like interference 
effect, a situation arises in which only the bonding or only the anti-bonding states is observed for 
a given photoemission geometry. 
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The assumption that electrons in deep core states do not participate in the bonding of 
sohds, and thus do not show band-hke dispersion, is at the base of our understanding of 
the sohd state. It is also of essential practical significance for many experiments, such as 
x-ray photoelectron spectroscopy. This technique derives its power from the fact that the 
precise value of the core level binding energy depends on the chemical environment of the 
emitting atoms, but it is tacitly assumed that it has a single, well-defined energy, i.e. it 
does not show any dispersion. Violations of this assumption have only been found for small 
molecules in the gas phase such as C2H2 or N2 with much stronger bonding and shorter 
bonding distances than present in solids [l, 2|. In this paper we report the observation of 
a considerable band-like dispersion of the C Is core level in graphene. The dispersion is 
observed as an emission-angle dependent binding energy modulation and it is shown that 
under appropriate conditions only the bonding or anti-bonding states can be observed. 

The binding energy modulations of the C Is core state are illustrated in Fig. [T]^a)-(e). 
Each panel shows a group of spectra taken at a fixed polar emission angle ^ as a function 
of azimuthal emission angle (f). Clear shifts of the peak position are observed. Fig. [H^f) 
shows a comparison of the spectrum taken at normal emission and one taken at 6' = 25° 
together with the result of a peak fit to these two spectra. Details of the fitting procedure 
are described in the supporting material. 

The binding energy variation obtained from the peak fitting is given in Fig. [T]^g) with 
the markers corresponding to the spectra in (a)-(e). Strong changes are evident with the 
largest difference of binding energies spanning ?«60 meV. The variation is consistent with 
the point symmetry of the graphene lattice. Fig. [Hh) shows the intensity variations of the 
peaks, displayed as the modulation function (see supplementary material). Strong inten- 
sity modulations are observed, caused by photoelectron diffraction in the final state. The 
variations follow the point symmetry of the graphene lattice, but they do not appear to be 
correlated with the binding energy in (g), with the main structures being in phase for some 
polar emission angles and out of phase for others. 

While the data of Fig. [1] serve to illustrate the nature of the effect and the fitting 
procedure, they are insufficient to pin down the physical origin of the modulation. Fig. [2] 
therefore shows a much more extensive data set measured over many polar and azimuthal 
angles and at different photon energies. Again, a good fit to all the spectra in the data set 
was obtained for a single C Is component using always the same line shape. Note that this 
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FIG. 1: (a)-(e) C Is photoemission spectra taken at a photon energy of 400 eV, for fixed polar 
emission angles in each panel but at different azimuthal emission angles (p (see sketch of the 
experimental geometry at the top of the figure) . The spectra are all normalized to the same height 
and shown as a group plot, such that binding energy variations become evident. The first and last 
spectrum in a range of azimuthal angles (p is indicated by a thicker line, (f) Comparison of the 
spectrum taken at ^ = 0° and one taken at 25°. The lines are the fits through the data points 
using the lineshape parameters described in the supplementary material, (g) C Is binding energy 
required to obtain a good fit for the curves shown above (markers) as well as for the entire azimuthal 
range measured (lines). The green horizontal line marks the binding energy at normal emission. 
The binding energy uncertainty is smaller than 10 meV. Error bars are not shown, (h) Intensity 
variation of the C Is peak as a function of azimuthal angle. The curves are shifted vertically for 
clarity. 
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excludes the existence of unresolved components, because their intensities would modulate 
differently, changing the shape of the peak. The left panel of the figure shows the resulting 
intensity modulation function while the right panel gives the binding energy modulation. 
The modulation function is compared to the simulation for a flat, free-standing layer of 
graphene. The agreement between experiment and calculation is excellent. 

The binding energy modulation, on the other hand, is shown as a function of ky, the 
wave-vector component parallel to the surface, which is the only relevant wave-vector for 
a two-dimensional system like graphene. As the portion of the reciprocal space covered 
by the experiment increases at higher photon energies, a periodic pattern emerges which, 
however, does not coincide with the reciprocal lattice mesh. At the origin (ky = (0,0)) the 
binding energy is always close to its maximum value. The experiments at hu = 400 eV 
and 500 eV show that the binding energy takes its maximum value also at the next-nearest 
neighbor reciprocal lattice points, while the experiments at 600 and 700 eV show that this 
occurs again at the next-nearest neighbors of these latter points. Interestingly, at all the 
other reciprocal lattice points the binding energy takes its minimum value. The periodic 
pattern, therefore, is described with two complementary sublattices: the binding energy is 
at its maximum value at all the points connected by vectors that are y/S times longer than 
the primitive vectors and rotated by 30°, and it is at its minimum value at all the other 
reciprocal lattice points. 

There are, in principle, several different mechanisms which could lead to the observed 
binding energy variations. The first is the existence of several unresolved components which 
change their relative intensity due to a final state effect (photoelectron diffraction) and 
thereby mimic a peak-shift. This appears highly unlikely in the present case, not only be- 
cause of the excellent agreement with the intensity calculation for a single component, but 
also because a peak shift of the observed magnitude would essentially require a complete 
supression of peak intensities in certain directions, and this is unrealistic for a usual pho- 
toelectron diffraction-type modulation. It would also be very unlikely that the mechanism 
leads to a binding energy modulation which is periodic in reciprocal space, as seen in Fig. |2l 

The second possible explanation is a recoil effect in which some of the photoelectron's 
energy is used to excite lattice vibrations. This effect has been reported for photoemission 
from the graphite C Is state [s^. For increasing photon energies it leads to a decreasing 
apparent binding energy since some of the photoelectron's energy remains with the emitting 
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FIG. 2: Left panel: stereographic projection of the photoemission intensity modulation as a func- 
tion of emission angle for scans taken at different photon energies hu. The colored fraction of the 
disk are the data, the greyscale part is a calculation of the expected intensity. Right panel: The 
corresponding binding energy variations obtained from the peak fitting. Note that these variations 
are shown as a function of the wave vector parallel to the surface rather than emission angle. The 
binding energy uncertainty is of the order of 10 meV. The green crosses correspond to the reciprocal 
lattice points of graphene, i.e. to the f points. 
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atom. The effect can be of considerable size, several hundred meV for photon energies of 
several keV, but its magnitude should be insignificant for the low energies used here. One 
would not expect it to lead to any periodic modulation either. 

This leaves a third possibility which is an initial state effect, i. e. a band-like dispersion 
of the initial state. The simplest conceivable picture for this is the formation of a a-type 
band between the Is states of the two atoms in the unit cell of graphene, highlighted in 
orange in the sketch of Fig. [1] A tight-binding calculation of such a band is shown in Fig. 
[3]^a). The absolute binding energy and the band width are arbitrarily chosen to mimic those 
observed here. The dispersion shows two bands with the highest energy separation at T and 
degeneracy at the K point of the two-dimensional Brillouin zone. In the following we refer 
to these bands, somewhat loosely, as the bonding and the anti-bonding band. 

On the face of it, it appears hard to reconcile such a dispersion with the experimental 
observation, as one would expect to observe a single, narrow C Is peak at K and a broad 
or even split peak at T. This is clearly not supported by the data. Also, the a-band should 
be periodic in reciprocal space, e. g. the peak position and width should be the same at all 
r points. According to the right panel of Fig. Ej this is not the case either. Most of the F 
points marked in the figure appear close to either a maximum or a minimum in the binding 
energy but clearly the observed periodicity is not the same as that of the reciprocal lattice. 

The hypothesis of band dispersion and the experimental data can, however, be reconciled 
when taking into account a curious interference effect which is caused by the presence of 
two atoms in the unit cell of graphene. In the most simple picture, this interference can be 
interpreted as a type of Young's double slit phenomenon in which the two atoms in the unit 
cell act as electron sources. The effect has been studied in detail for the valence band of 
graphite and graphene 0, 0] . 

As an illustration of the interference effect, we have calculated the expected photoemission 
intensities for all the anti-bonding and all the bonding states. The results are given in Fig. 
[3]^b) and (c), respectively. The strength of the interference effect is evident: in the first 
Brillouin zone, for instance, emission from the anti-bonding states is entirely suppressed 
while it is intense from the bonding states. In the neighboring zones it is the other way 
round. Figs. |3](d)-(f) show the emission intensity for smaller energy windows at the top of 
the anti-bonding band, at the bottom of the bonding band, and just below the Dirac point 
at K. In the last case, a constant energy contour shows a triangular shape, as expected for 
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FIG. 3: (a) Result of a tight-binding calculation for a cj-type band formed from the C Is core states 
in graphene. The bonding (blue) and anti-bonding (red) bands are degenerate at K and show the 
highest splitting at f. The inset shows the Brillouin zone of graphene. (b) and (c) Calculated 
photoemission intensity from all the anti-bonding and bonding states, respectively. The grey-scale 
is chosen such that bright corresponds to high intensity. The green crosses mark the reciprocal 
lattice of graphene and the green hexagon the first Brillouin zone, (d)-(f) Calculated photoemission 
intensity from the states in the binding energy windows indicated by the small circles in (a). Note 
the similarity of the emission pattern from the bonding states in (d) with the positions of maximum 
binding energy in the right panel of Fig. [2j 

the non-linear dispersion away from K. The photoemission intensity around this triangular 
contour shows strong variations which are caused by the interference effect and very similar 
to the results obtained for the valence vr-band of graphite and graphene 0, 0]. 

As already expected from Fig. [3t^b) and (c), the interference effect is even stronger for 
emission from the bonding and anti-bonding states at F (see Fig. EJ^d) and (e)). The 
intensity variations of both bands are almost opposite to each other: for some F points only 
the bonding band is observed, for others only the anti-bonding band. For normal emission 
this is easy to understand: for the bonding band the wave functions centred on the two atoms 
in the unit cell emit in phase and this band is observed. For the anti-bonding wave function, 
the two atomic wave functions emit out of phase, thus suppressing the photoemission. 
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The presence of the interference effect easily reconciles the hypothesis of a cr-band for- 
mation with the data. First of all, it explains the fact that the peak shape is similar for all 
emission directions. For emission near K the peak is narrow because of the degeneracy of 
the bands at K . At the F points the peak is also narrow, in contrast to naive expectation, 
because it does not show both the bonding and the anti-bonding bands but rather only one 
of them at every given F point. The interference effect is also responsible for the apparently 
incorrect periodicity of the band dispersion in reciprocal space. The situation is not such 
that all F points are equivalent, because one either observes emission from the bonding or 
from the anti-bonding band. Note that the sign of the observed modulation is consistent 
with this interpretation, too. The binding energy of the peak observed at normal emission 
is close to the global maximum of the entire data set, consistent with emission from the 
bonding cr-band. In fact, the calculated emission from the bonding band shown in Fig. [3]J^d) 
can be compared directly to the plot of the apparent binding energy in Fig. [21 which is also 
scaled such that high binding energies are bright. 

The size of the bonding / anti-bonding splitting in graphene can be inferred from the 
difference of observed binding energies at inequivalent F points. Combing all the available 
data from F points showing either the bonding or the anti-bonding band, we evaluate the 
size of the splitting to be 60 ± 10 meV. We can compare this value to the size of the bonding 
/ anti-bonding splitting in small carbon-containing molecules such as C2H2 (C-C distance 
of 1.2 A, splitting of 105 meV |l^) and C2H4 (C-C distance of 1.34 A expected splitting of 
20-30 meV jsl) if we assume that the matrix element for hopping between the core electrons 
on the two carbon atoms depends exponentially on the C-C distance and that the size of 
the splitting scales linearly with the number of nearest neighbors. Assuming a splitting of 
105 meV and 25 meV for C2H2 and C2H4, respectively, we would expect a total bandwidth 
of ~30 meV for graphene. The observed magnitude of the dispersion is somewhat larger but 
in the right order of magnitude. 

Apart from the fundamental importance of our results for the bonding in solids, a de- 
pendence of the core level binding energy on the emission angle can have implications for 
the interpretation of high-resolution core level data from graphene, graphite and related 
materials. The absolute magnitude of the binding energy variation observed here is appre- 
ciable compared with usual chemical shifts and could easily be interpreted incorrectly. In 
a limited data set, for instance, an apparent shift of the C Is peak as a function of polar 
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emission angle might be mistaken for signs of a surface core level shift in graphite because of 
the higher surface sensitivity for off-normal emission. Ignoring dispersion effects might play 
some role in the recent dispute about the existence of a surface core level shift in graphite 
7H10| , but note that the interference effect reported here is inconsistent with the observation 
of multiple C Is components. Here we can state that to a good approximation either the 
bonding or the anti-bonding band can be observed in any given emission direction. 

In conclusion, we have shown that the interaction between the two atoms in the unit 
cell of graphene is sufficiently strong to induce the formation of a a-band derived from 
the carbon Is state. Interference effects in the photoemission process give rise to a strong 
suppression of either the bonding or the anti-bonding band for different emission directions, 
such that the absolute size of the splitting could be determined from the entire data set. 
The interference effect furthermore allows to observe only the bonding or the anti-bonding 
states under certain conditions, opening the unique opportunity for detailed and possibly 
local as well as time-resolved studies of bonding in sohds. 
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SUPPORTING MATERIAL: EXPERIMENTAL AND CALCULATIONAL METH- 
ODS 



The graphene film on Ir(lll) was prepared in situ using standard procedures [HI] and 
the sample quality and cleanliness was monitored by low energy electron diffraction (LEED) 
and x-ray photoelectron spectroscopy. Angle-resolved photoemission data were taken at the 
SuperESCA beamline of the storage ring ELETTRA. The experimental chamber is equipped 
with a 150 mm electron energy analyser from SPECS, implemented with a delay line detector 
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developed in-house. It allows fast acquisition of the spectra in snap-shot mode, with up to 
thousand data points in each spectrum and in as short as 80 ms/spectrum. The data shown 
in Fig. m are acquired in scanning mode at hu =350 eV, 400 eV and 700 eV, while in 
snap-shot mode at 500 eV and 600 eV. The fast scans, with a total of 420 spectra that fill 
1/3 of the emission hemisphere, take less than 30 min. At the other photon energies more 
spectra were taken: 2270 at hu = 700 eV, 820 at hu = 400 eV and 1300 at hu = 350 eV. 
Such scans take several hours. The photon energy was calibrated by the absorption edges 
of gas-phase Ar, Ne and N2. The manipulator used in the experiments is a modified version 
of the VG-CTPO with fully computer-controlled polar and azimuthal rotations. 

For the analysis of the data all the spectra taken at one photon energy were first aligned to 
the Fermi level of the Ir substrate, then they were fitted with the same line shape parameters 
but the binding energy was left free in the fit. Fits for two spectra are shown in Fig. 
The actual peak fits were made using a Doniach-Sunjic line profile [l^ with a Lorentzian 
width of 130 meV, an asymmetry parameter of 0.093, a Gaussian width depending of the 
photon energy and a linear background. The narrowest Gaussian width was 165 meV and 
this could be used for photon energies up to 400 eV. It is the value used for the two fits 
shown in Fig. [H Note, however, that the spectra were measured and fitted in a wider range 
than displayed in Fig. [1] We have observed a small variation of the lineshape over the entire 
data set taken at any given photon energy. However, a good fit to all the spectra can be 
obtained with one set of parameters describing the line shape. Arbitrariness in the exact 
values of these parameters induces a small uncertainty in the absolute degree of binding 
energy modulation of about 5 meV. We thus estimate the total uncertainty of the observed 
bandwidth to be of the order of 10 meV. 

The modulation functions in Fig. [T](h) and the left part of Fig. [2] were obtained for 
each polar emission angle 9 from the peak intensity J(9,0) as (1(6,0) — /o(6))//o(6), 
were /o(0) is the average value of each azimuthal scan. Note that no artificial symmetry 
was imposed on the data. The experimental modulations are as symmetric as they appear, 
greatly enhancing out confidence in the data. The calculated modulations in Fig. |2] were 
obtained in the same way. 

The potoemission intensity calculations used to obtain the calculated modulation func- 
tions in Fig. |2] were performed by the program package for Electron Diffraction in Atomic 
Clusters (EDAC) [l3|. The photoemission intensity at each emission angle was calculated 
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as the incoherent sum of the intensities for the two atoms in the graphene unit cell. This 
approach might, at first sight, appear inconsistent with the main conclusion of this paper 
that the Is state is described as an extended coherent wave function. Note, however, that 
what we calculate and compare to is the sum of the bonding and anti-bonding states. The 
sum of the probability density for the wave functions is strongly peaked at the atomic cores, 
independent of the wave vector chosen. Consequently, one would expect the same photoe- 
mission signal, including final-state diffraction effects, as from the sum of all core states. 
The simulations were performed on a flat, free-standing graphene layer. The lattice param- 
eter was set to 2.466 A. The influence of the underlying Ir(lll) substrate, that of small 
changes of the lattice parameter, as well as that of a relatively small corrugation of the 
graphene layer with moire periodicity [14] was tested. No significant changes were found in 
the diffraction patterns. Parameters for the calculations are: Vq = 8 eV, linear polarization, 
angle between incoming light and analyzer 70°, angular acceptance 10°, system at room 
temperature. The Debye temperature of graphene is much higher than this and hence its 
precise value does not lead to any noticeable differences in the result. As seen in Fig. [2] 
the agreement between these calculation and the data is excellent with the two being nearly 
indistinguishable. However, in order to obtain this agreement, the photoelectron's kinetic 
energy in the simulation had to be assumed to be 12 eV lower than that in the experiment. 
The reason for this is presently unknown. 

The tight-binding calculations for the band structure and photoemission intensity in 
Fig. [3] are very similar to those presented for the vr- valence band in Ref . js | , but use cr-band 
symmetry instead. The parameters used for the calculation are inspired by the experimental 
binding energy and band width but they are only chosen for illustrative purposes. Their 
precise value is irrelevant for the conclusions of this paper. 
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